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Algorithms for quantum information processing are usually decomposed into 
sequences of quantum gate operations, most often realized with single- and two- 
qubit gates [lj. While such operations constitute a universal set for quantum 
computation, gates acting on more than two qubits can simplify the implemen- 
tation of complex quantum algorithms [2]. Thus, a single three-qubit operation 
can replace a complex sequence of two-qubit gates, which in turn promises faster 
execution with potentially higher fidelity. One important three-qubit operation 
is the quantum Toffoli gate which performs a NOT operation on a target qubit 
depending on the state of two control qubits. Here we present the first experi- 
mental realization of the quantum Toffoli gate in an ion trap quantum computer. 
Our implementation is particular efficient as we directly encode the relevant logic 
information in the motion of the ion string. 

Quantum gate operations have been realized with several physical systems, e.g. nuclear 
magnetic resonance (NMR)|i3|, single photons[U E], superconducting devicesjHl El E], 
cavity QEDfTDj. or ion traps [HI [T21 EH EH]- All quantum gates demonstrated so far are 
two-qubit entangling operations, apart from one implementation of a three qubit Toffoli 
gate in NMR[15]. For ion traps several schemes for implementing a quantum Toffoli gate 
have been proposed[16j [17]. The quantum Toffoli gate is valuable in complex quantum 
algorithms and has an immediate practical application as correcting operation in quantum 
error correction [151 HH]- We report here on a realization of the Toffoli gate with an ion trap 
quantum computer which requires half the resources compared to implementations based on 
concatenated two-qubit gates [19]. This significant contraction is achieved by encoding the 
combined information of the two control qubits in the common vibrational mode of the ion 
string. 

Our experimental system consists of a string of 40 Ca + ions confined in a linear Paul trap. 
Each ion represents a qubit, where quantum information is stored in superpositions of the 
Si/2(m=— 1/2) = \S) = |1) ground state and the metastable D 5 / 2 (m=— 1/2) = \D) = |0) 
state of the 40 Ca + ions[20j. The centre-of-mass (COM) vibrational mode of the ion string 
is used to mediate the interaction between the ion qubits. Each experiment includes a) the 
initialization of the qubits and the COM mode in a well defined state, b) the actual gate 
operation and c) a quantum state measurement. Initialization of the COM mode to the 
ground state is achieved by Doppler cooling on the Si/ 2 <-»■ P1/2 dipole transition followed 



by sideband cooling on the qubit transition. Optical pumping ensures that all qubits are 
prepared in Si/ 2 (m=— 1/2). The gate operation consists of a series of laser pulses which 
are applied to individually addressed ions. The electronic and vibrational state of the ion 
string are manipulated by setting frequency, length, intensity, and phase of the pulses. The 
final state of the ion qubits is measured by scattering light at 397 nm on the S1/2 <-► P1/2 
transition and detecting the fluorescence with a CCD camera. The presence or absence 
of fluorescence of an ion corresponds to a projective measurement in the {\S), \D)} basis. 
Further details about the experimental setup can be found in Ref. [2"U] . 

The pulse sequence for our implementation of the quantum Toffoli gate consists of three 
major parts: 1) Encoding of the joint quantum information of the control qubits \c\) and 
|c 2 ) in the vibrational COM mode, 2) a NOT operation on the target qubit controlled by 
the vibrational mode, and 3) the reversal of the encoding step 1). Here, the encoding is the 
crucial part of the gate. We exploit the fact that it is always possible to add more phonons 
to the vibrational mode, while phonons can only be subtracted until the vibrational ground 
state is reached. The joint encoding of the logic control information (ci AND c 2 ) into the 
vibrational mode can be depicted as a two-step process: Initially, the vibrational mode 
contains no phonons, i.e. \vib) — \n — 0). In the first step, two phonons are deposited in 
the COM mode if both control ions are in \c1C2) = \SS) while less than two phonons are 
generated for the other three basis states (see Methods). In the second step, one phonon 
is removed such that only for the initial state \c\C2) = \SS) one phonon remains in the 
COM mode. Thus, we encode the logic information (ci AND c 2 ) in the vibrational mode as 
(ci AND c 2 ) = 1 =>• \n — 1) and (c\ AND c 2 ) = =>- \n = 0). The implementation of the 
second step requires the use of a composite pulse sequence which realizes the subtraction of 
one phonon regardless of the initial phonon number (for details see Tab. II] and Methods). 
The next part of the Toffoli gate consists of a NOT operation on the target qubit depending 
on the binary information encoded in the vibrational state. This operation is provided by the 
central part of the Cirac-Zoller controlled NOT gate[2T] as demonstrated in[13l[Tlj. Finally, 
the third part reverses the encoding and restores the vibrational mode to the ground state. 



The ideal unitary map implemented by our pulse sequence, see Table [TJ is given by 
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The factor exp(—in jj^a z,t) , with the Pauli operator az, is a phase shift on the target qubit 
\t). This phase factor can either be included in a redefinition of the computational basis or 
compensated for by single qubit operations. In total, this unitary evolution is fully equivalent 
to an ideal Toffoli gate. 

The unitary map L^p describes an operation which applies a NOT operation on the target 
ion if and only if the control ions are in the state \SS) . We verify this behaviour by preparing 
the system in all eight basis states and by measuring the output state probabilities after the 
Toffoli operation. Figure [I] shows a bar graph of these probabilities, the so-called truth table. 
Ideally, this shows the absolute values of the respective elements of the unitary matrix. On 
average we obtain probabilities of ~ 81% that the ion ends up in the correct output state. 
However, this analysis only reveals the classical features of the gate. In order to prove that 
we indeed perform the desired quantum operation, a more thorough analysis is required. 

A complete characterization of the action of our pulse sequence must take into ac- 
count experimental imperfections, in particular the coupling of our qubits to their envi- 
ronment, which will in general lead to a non- unitary time evolution. Such a characteriza- 
tion is provided by the operator sum representation in which a process £ is described by 
Pout = £{pin) = Em,n=i Xm^AnPinA^ , where p in and p out are the input and output density 
matrices, iV is the number of qubits, and A is a basis of operators in the Hilbert space of 
dimension 2 N x 2 N [22]. The process matrix x contains all information about the inves- 
tigated quantum operation and is obtained by quantum process tomography (HJ [23]. The 
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FIG. 1: Experimentally obtained truth table of the Toffoli gate: After preparing the three ions in 
one of the eight basis input states \SSS) to \DDD) (right axis) the probabilities of all basis output 
states (left axis) are measured 100 times. On average the correct output state is reached with a 
probability of 81(±5)%. The error mainly originates from quantum projection noise. 



absolute values of both the ideal and the experimentally obtained x-matrices are depicted 
in Fig. [2} The "fingerprint" of the ideal x-matrix is clearly visible in the experimental data. 
A quantitative measure for the performance of the gate operation is the mean gate fidelity 
Fmean = mean^,. (ipi\Urp £(\ipi)(ipi\) f/'jlV^) , where £Ap is the ideal unitary map and the 
ipi are a set of 2 x 10 5 pure input states randomly drawn from the Haar measure for the 
unitary group £/(8)|24j. The mean fidelity is the average value of the fidelity between the 
output state expected from the ideal unitary transformation and the output state calculated 
from the experimentally determined x matrix. For the results shown in Fig. |2j we ob- 
tain F mean = 71(3)%, where the error is due to statistical uncertainties in the tomographic 



measurements [25] . The observed infidelity is mainly caused by technical imperfections. Ad- 
dressing errors due to the finite width of the focused laser beam result in a weak excitation 
of the neighbouring ions. For the given setup, the addressing error described by the ratio 
of non- addressed to addressed Rabi frequency was on the order of 7%, causing an infidelity 
of about 12%. Temperature changes and voltage fluctuations of the trap electrodes produce 
a drift of the vibrational mode frequency. Assuming a detuning of 100 Hz, this leads to an 
infidelity of 7%. Note that these two errors are purely technical imperfections and do not 
represent physical limitations. Remaining infidelity contributions are: initialization of the 
COM mode in the ground state - about 1%; laser linewidth and magnetic field fluctuations: 
about 1 %; ion state initialization: about 0.5% per ion. The decoherence time of quantum 
information stored in the vibrational COM mode is on the order of 85ms, while the heating 
rate is on the order of 1 phonon per 140ms. These two effects do not limit the performance 
of the quantum Toffoli gate at the current stage. 

A comparison of this result to the fidelity expected from a realization based on two- 
qubit CNOT gates highlights the benefits of our scheme. With our experimental setup, 
two-qubit CNOT gates have been demonstrated [2] with a mean fidelity of F mean = 92.6%. 
We estimate that a realization of a composite quantum Toffoli gate with six such CNOT 
gates [23] would result in a gate fidelity of approximately F = (92. 6%) 6 ~ 63%. We note 
that this is an optimistic estimate as it relies on two-ion data only. The duration of our gate 
is 1.5 ms, which is mainly determined by the coupling strength for transitions of the kind 
\S,n) ^ \D,n + 1) (see Methods). A two-qubit controlled NOT gate takes approximately 
700/is [Hj. Therefore, realizing a Toffoli gate based on six CNOT gates [23] would be about 
three times longer than the sequence described above. Finally, an adaption of the proposal 
for the Toffoli gate using qudits [T7] to an ion trap system would also be shorter than a 
decomposition into CNOT gates, but still 30% longer than our sequence. 

To conclude, we have shown the first implementation of a three-qubit Toffoli gate with 
trapped ions. It is notably faster than any other proposed scheme. The achieved fidelity 
of F mean = 71(3)% is higher than an optimistic estimate of the Toffoli gate fidelity based 
on two-qubit CNOT gates. The particular efficiency of this scheme is achieved by directly 
encoding the relevant logic information of the control qubits in the harmonic oscillator states 
of the vibrational mode. Such encoding may also help to simplify other complex quantum 
operations and may be applied to any system in which qubits are coupled via a harmonic 



oscillator. In addition, the Toffoli gate has an immediate practical application in quantum 
error correction. Here, depending on the state of the qubits carrying the error syndrome a 
conditional bit flip has to be applied to the qubits carrying the information [To! [T5] . This task 
can either be accomplished by a projective measurement followed by classically controlled 
operations or directly by the Toffoli gate. The latter approach is particularly favourable in 
quantum systems where a projective measurement is not feasible or not efficient. 

Methods 

Single qubit operations are implemented by laser pulses at frequency u> resonant with 
the transition \S,n) <-> \D, n), where n is the phonon number of the state. The motion of the 
ion string is manipulated by laser pulses at frequency ujq + u z , where uj z is the COM mode 
frequency along the trap axis, inducing transitions of the kind \S,n) <-> \D,n + 1). This 
allows us to remove a phonon from the system, provided an ion is in the \D) state. Note that 
the Rabi frequency on transition \S,n) <-> \D, n + 1) is Q n = Q n=0 \/n + 1 (Q = 2tt 3.3kHz). 
During the encoding section, the first two pulses of the Toffoli pulse sequence (see Table [I]) 
map the state of the control qubits on the motion: 

\SS,0) -► \DD,2) (1) 

\DS,0) - sin *\DD,1) + cos ^\DS,0) (2) 

\SD,0) -> cos^=\DD,l)-sm^=\DS,0) (3) 

\DD,0) -> 1^0,0,0) (4) 

Then the first control bit is always in state \D). The next step is to remove exactly one 
phonon by transferring control ion 1 from \D,n + l) to \S,n). This is not directly possible in 
all four cases due to the different coupling strengths on the \D, 2) <^ \S, 1) and \D, 1) <-> \S, 0) 
transitions. However, a sequence of pulses of different length and phase can be tailored such 
that in spite of the different coupling strengths the operation \D,n+ 1) <-► \S, n) for n — 1, 2 
is realized. The mechanism of this composite pulse sequence - pulses 3-5 in Table [I]- can be 
understood using the Bloch sphere picture with state \D, n + 1) as the north pole and IS*, n) 
as the south pole (e.g. see Ref. |20J). Starting from \D, 2), pulse 3 flips the state vector from 
the north pole to the equator, pulse 4 rotates the vector around itself, and pulse 5 continues 
the flip to the south pole reaching 15,1). On the other hand, for the initial state \D, 1), 
pulse 3 rotates the state vector only partially towards the equator, pulse 4 mirrors the state 
vector on the equatorial plane, so that pulse 5 can complete the rotation to the south pole 



\S, 0). For zero phonons nothing changes and thus conserves \D,0). After this step a single 
phonon is remaining only for the initial control state \SS). 
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FIG. 2: Comparison between (a) ideal and (b) experimental process x-matrices of the Toffoli 
gate. The absolute values of the elements are shown in the operator basis a c \ (g) a C 2 <S> at £ 
{I ® I ® 1,1 <8> I <g> X, . . . ,Z <& Z <& Z}. Ideal and experimental process matrices show the same 
"fingerprint" . 
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Encode first target qubit onto motion 


Encoding 
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D+ ( 7T 37T \ 
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Encode second target qubit onto motion 


3 

4 
5 


i?+(7T,0) 


Composite pulse to remove one phonon 
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Prepare target qubit for motion controlled NOT 


controlled NOT 


7 
8 
9 


i? 3 + (V2n, f ) 

^(f,o) 


Composite phase gate 


10 


«b(|,(^-1K) 


Complete motion controlled NOT on target qubit 


11 
12 
13 


Rt (*, (-1 + ^K) 


Undo encoding algorithm 


Decoding 


14 


^fed + ^k) 


Decoding finished for second control qubit 


15 


^M + ^V) 


Decoding finished for first control qubit, Toffoli complete 



TABLE I: Pulse sequence of the Toffoli gate implementation: Laser pulses on the ith ion on the 
15, n) <-» \D, n) transition are denoted by Ri(9, <ft), on the IS", n) <-> \D, n + 1) transition by Rf{6, </>), 
with laser phase (j) and pulse area = Qt in terms of Rabi frequency $7, and pulse length t, for 
details see ref. |20] . 
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